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We investigate the Josephson effect in one triple-terminal junction with embedded parallel-coupled 
double quantum dots. It is found that the inter-superconductor supercurrent has opportunities to 
oscillate in Ai r period, with the adjustment of the phase differences among the superconductors. 
What is notable is that such a result is robust and independent of fermion parities, intradot Coulomb 
strength, and the dot-superconductor coupling manner. By introducing the concept of spinful 
many-particle Majorana modes, we present the analytical definition of the Majorana operator via 
superposing electron and hole operators. It can be believed that this work provide a simple but 
feasible proposal for the realization of Majorana modes in a nonmagnetic system. 

PACS numbers: 71.10.Pm, 73.23.-b 


Majorana fermions, exotic quasiparticles with 
non-Abelian statistics, have attracted enormous at¬ 
tention because of their fundamental physics and 
potential application in topological quantum com¬ 
putation. During the past years, various schemes 
have been proposed to realize unpaired Majorana 
fermions, such as in a vortex core in a p-wave super¬ 
conductor [1-5] or superfluid [6-8]. Feasible propos¬ 
als are to fabricate Majorana bound states (MBSs) 
at the ends of a one-dimensional p-wave supercon¬ 
ductor, which has been constructed by adhering a 
Rashba nanowire with a strong magnetic field or a 
ferromagnetic chain to the s-wave superconductor 
[9-14]. However, how to verify the existence of MBSs 
is a key issue and is rather difficult. 



FIG. 1: (a) Schematic of a Josephson junction with a 
DQD molecule, (b) Difference between the GS levels in 
even and odd FPs. (c)-(d) Josephson currents in the 
even and odd FPs, respectively. 


Superconductor-quantum-dot systems are also 
proposed to realize Majorana modes, by including 
spin-orbit coupling or anisotropic magnetic fields 
[15-17]. Although quantum dots (QDs) have several 
advantages, strong spin-orbit or anisotropic mag¬ 
netic fields form a major hurdle and limit ma¬ 
terial flexibility. For example, anisotropic mag¬ 
netic fields only result in spinless localized Ma¬ 


jorana modes when Ez ^ A ~ t with A the 
superconducting gap and t the intersite hopping 
[17]. In order to overcome such difficulty, the struc¬ 
ture of superconductor/double-QD/superconductor 
(SC/DQD/SC) has recently been designed [18]. It 
has been observed that in such a structure, spinful 
many-particle Majorana modes can arise in the pres¬ 
ence of appropriate Coulomb interactions and local 
magnetic fields applied on the QDs. However, for 
practical experiments, it is a formidable challenge 
to apply two different local magnetic fields on the 
nanoscale. Therefore, any new suggestion to realize 
the Majorana modes in the non-magnetic systems is 
desirable. 


Some recent researches suggest that a triple¬ 
terminal Josephson junction can be a promising can¬ 
didate for achieving Majorana modes. Main reason 
is that in this junction, the Kramers degeneracy of 
its different-parity energy levels can be broken in the 
presence of superconducting phase difference, and 
adjusting phase differences can change the ground- 
state (GS) parity [19]. Motivated by these phenom¬ 
ena, in this Letter we would like to investigate the 
Josephson effect in one triple-terminal junction with 
embedded parallel-coupled DQDs. As a result, by 
adjusting the phase differences among the supercon¬ 
ductors, the 47r-periodic oscillation of the inter-SC 
supercurrent appears, which is independent of struc¬ 
tural parameters and fermion parities. Then, we de¬ 
fine the Majorana operator via superposing electron 
and hole operators. Hence, this work provides a sim¬ 
ple but feasible proposal to realize Majorana modes 
in non-magnetic systems. 


The considered Josephson junction is illustrated 
in Fig. 1(a), where parallel-coupled DQDs couple to 
three SCs, respectively. The Hamiltonian for the 
electron motion in this system can be written as H — 
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FIG. 2: (a) Difference between the GS levels in even and 
odd FPs in the case of infinite Coulomb interaction, (b) 
Corresponding Josephson current in the even-FP case. 
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a a/ccr ( a aka) is an operator to create (annihilate) an 
electron with momentum k and spin orientation cr in 
SC-<x denote the corresponding energy levels. 

A is the superconducting gap and is the phase of 
the superconducting order parameter. d}- a ( dj a ) is 
an operator to create (annihilate) an electron with 
energy Ej and spin orientation a in QD-j, respec¬ 
tively. Uj indicates the strength of intradot Coulomb 
repulsion, and t c is the interdot coupling coefficient. 
V a u(D ) denotes the coupling between SC-<a and QD- 
1 ( 2 ). 

In this work we would like to take a few simplifica¬ 
tions to make clearer the Josephson effect in such a 
SC/DQD/SC structure. The SCs are assumed to be 
identical (e a k = £k and A a = A) except for a finite 
phase difference ipt = —(fR = cp/2 and <pm = 0. The 
Josephson current between at zero temperature can 
be evaluated by deriving the GS energy Eqs with 
respect to the superconducting phase difference, i.e., 

2ed E GS (ip) 

= (2) 

It should be emphasized that since we are inter¬ 
ested in transport processes involving Andreev re¬ 
flection, we focus on the limit A —>> oo. In this 


case, as the quasiparticles in the SCs are inacces¬ 
sible, tracing out of the degrees of freedom of the 
SCs does not induce any dissipative dynamics in the 
DQD system and it can be performed exactly. The 
resulting Hamiltonian dynamics of the DQDs can 
be described by the following effective Hamiltonian 
[ 20 ] 
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^aU(D), defined by T aU{D) = \V aU{D) \ 2 5(u - 

Sk), describes the coupling strength between the 
DQDs and leads. Within such an approximation, 
the Josephson current between two SCs can be eval¬ 
uated by diagonalizing It is known that in 

this structure, fermion parity (FP) is good quantum 
number. The Hamiltonian matrix should be written 
out according to FPs, namely, 
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with T af3 = To and A = To(l + 2 cos ^). Following 
the theory in the above section, we proceed to dis¬ 
cuss the property of the Josephson current between 
SC-L and SC -R. 

In Fig. 1(b) we define E g = Eq^—Eq^ and present 
the difference between the different-FP GS energies. 
The structural parameters are taken to be ej = £, 
t c = 4, and U = 10, respectively. The result in this 
figure shows that in such a structure, the Kramers 
degeneracy of the GS energies in different GSs is 
clearly broken. Changing e and ip can alter the GS 
parity, with the transition position (i.e., position of 
Eg = 0) denoted by the dotted lines. In addition, 
one can find that in most of the parameter space, 
the energy of the odd-parity GS is lower than the 
even-parity GS. Such a result is opposite to that in 
the normal SC-existed system where E g is always 
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less than zero due to the presence of Cooper pairs 
[19]. The other phenomenon is that E g presents dif¬ 
ferent oscillation manners in the regions e > 0 and 
5 < 0. Following the result in Fig. 1(b), we plot 
the spectra of the Josephson currents in different 
FPs, as shown in Fig.l(c)-(d). It is evident that 
both of the Josephson currents in different FPs os¬ 
cillate in 47r period with the similar directions of 
them. However, the differences between I e and I Q 
are easy to observe. In addition to their differ¬ 
ent magnitudes, the odd-FP current is suppressed 
completely in the region 5 < —7.0. This result 
can be explained by solving equation = 

i^( 0 ) |J/(°)). Under the condition of the parameters 
in Fig.l, the eigenvalues of can be listed, i.e., 
4°i) = 6 + 36, £$) = 4 + g, and = £ ff a) = 

5 + 2e =b ^93 + 185 + 5 2 + 16 cos f + 8 cos p, respec¬ 
tively. Note that the double degeneracy of these 
eigenvalues originates from the spin symmetry of ba¬ 
sis of . Surely, with the decrease of 6, E^ have 
an opportunity to become the GS level. And then, 
its independence of p suppresses the Josephson cur¬ 
rent directly. As is known, the degeneracy of the 
different-FP GSs allows to define a pair of Majorana 
operators (71, 72) which transform the two GSs into 
each other. Therefore, we ascertain that the Majo¬ 
rana modes can be realized in such a triple-terminal 
junction. 


Coulomb interaction is a key factor to influence 
the properties of Josephson current. As an extreme 
case, when V increases to infinity, the double oc¬ 
cupation in each QD will be forbidden, which will 
inevitably transform the Josephson effect. Next, we 
would like to discuss the Josephson effect by tak¬ 
ing U —)► oc. In such a case, the dimension of 
will decrease to five, and its eigenvalues are E = 
i?2(3) = 25, and = 5 d= ^/6 + 5 2 + 8 cos | + 45, 

respectively. On the other hand, will reduce to 
one 4x4 matrix with the eigenvalues E^) = s + t c 

and Eg°^ = e — t c . With these results, we can 
write out the expression of E g directly, i.e., E g = 
— y/6 + 5 2 + 8 cos f + 4 cos p + t c . Evidently, E g is 
a even function of 5. Its detailed property can be ob¬ 
served in Fig.2(a). Suppose E g = 0, the transition 
position of the GS parity can be well-defined with 
the relationship 5 = d= ^/t 2 — 6 — 8 cos f — 4 cos p. 
As for the Josephson current in the limit U —)> 00, 
it comes into being only in the even-parity case be¬ 
cause of the p independence of E^°\ As shown in 
Fig.2(b), the Josephson current oscillates in 47r pe¬ 
riod and its maximum shifts to the point of 5 = 0, 
though the current magnitude is somewhat weak¬ 
ened in comparison with the case of U = 10. 


The results in Fig.2 show that Majorana modes 
survive in the extreme case of U —>• oc. Due to the 



FIG. 3: Even-FP Josephson current and GS level as 
functions of p in the asymmetric QD-superconductor 
coupling case, (a)-(b) Josephson current and GS level 
in the case of up-down asymmetry with Flu = F ru = 
To + S and Fld = Frd = To — S. (c)-(d) Results of 
left-right asymmetry with Flu — Fld — To + J and 
Fru = Frd = To — S. (e)-(f) Results of up-down 
left-right asymmetry with Flu = Frd = To + S and 
Flu — Frd = To — 8. 


decrease of the dimension of in such a case, we 
are able to present the definition of the Majorana 
modes in an analytical way. The role of Majorana 
operator is to change the GS parity, which can be 
expressed as 


7i|^(°)) = |^ (e) ). 


(4) 


Taking the case of p = tt as an example, where 

= V2+(U/2+i+ 2 [_(£ + v/ 2+ i? )l°°) - I 

> + ltl>] and = |[| t 0) - |0 f) + U 0) - |0 i)], 
the Majorana operator can be worked out as follows. 
First, we define the spinful Marjoana mode with its 
operator being 71 = an u 7^ + + cn 2 4,711 + 

dri2x 7 it + en 2t7u + f n 24714 + # 7 it + ^724. + jn u n 2 ^ 

and 7 ja- = -E(dt a + dj a ). Next, by using the re¬ 
lationship 71 = |\lA e ))(T( 0 )|, the relevant parame¬ 
ters in the expression of 71 can be obtained. They 

,- , , b = C = 

V2+G+v / 2+72)2 


are a = d = -g = 


2^/2(s-\-y/2-\-s 2 — 1) 


, and e = 


2 V 2 


whereas 


V2+( £ +v / 2+7^P ’ V /2+(e+72+7^P 

f = h = j = 0. It is certain that 72 has a similar 
form with site indexes 1 and 2 interchanged. Based 
on the above analysis, we ascertain that in such a 
triple-terminal DQD junction, the Majorana mode 
can be constructed. Since the realization of Majo¬ 
rana mode is irrelevant to the magnetic factors and 
other fields, such a scheme is more feasible compared 
with the previous works. 

Since the parallel geometry of such a Josephson 
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junction, the symmetry of QD-SC couplings can be 
adjusted to be various manners, such as up-down 
asymmetry, left-right asymmetry, and up-down left- 
right asymmetry. Different asymmetry manners will 
modify the Josephson effect in special ways. We 
next suppose V —>• oc and discuss the influence of 
the asymmetry of the QD-SC coupling manner on 
the formation of Majorana mode. The results of 
the Josephson current and GS level are shown in 
Fig.3. It can be found that regardless of the QD- 
SC coupling manners, the increase of 8 efficiently 
decreases the magnitude of the Josephson current, 
followed by its weakened oscillation near the posi¬ 
tion <p = 2tt. At the asymmetric limit of 8 = 0.9, 
the additional oscillation of I e around ip — 2n disap¬ 
pears, as a result, the Josephson current obeys the 
relationship that I e ~ sin ^ [See Fig.3(a), Fig.3(c), 
and Fig.3(e)]. Such a result can be well understood 
based on the GS-level oscillation shown in Fig.3(b), 
Fig.3(d), and Fig.3(f). These results mean that the 
case of asymmetric QD-SC coupling can not destroy 
the Majorana mode but promote its appearance. 

In conclusion, we have studied the Josephson ef¬ 
fect in one triple-terminal junction with embedded 
parallel-coupled DQDs. As a consequence, it has 
been found that the inter-SC supercurrent oscillate 
in 47r period, following the adjustment of the phase 
differences among the SCs. Moreover, such a result 
is robust and independent of FPs, intradot Coulomb 
strength, and the QD-SC coupling manner. By in¬ 
troducing the concept of spinful many-particle Ma¬ 
jorana modes, the analytical definition of the Ma¬ 
jorana operator has been presented. In addition, 
we have found that the asymmetric QD-SC coupling 
manner is conducive to the realization of Majorana 
modes in such a system. In view of all the results, 
we consider that this system can be a promising can¬ 
didate for the realization of Majorana modes in a 
nonmagnetic system. 
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